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Abstract 

Relays are used industrially in a wide variety of applications. 
Traditional mechanical relays are large, slow, noisy devices, 
but are still widely used in various machines and processes 
for control purposes. In this paper cross section of a moving 
electromagnetic relay is A; can be inside the external vertical 
hole magnetic chamber. Gap between the two is negligible, r 
is the resistance in ohms, and the coil has N turns. Once a 
potential source v volt DC and again, AC power is applied to 
the ends of the coil. Magnetic material can be assumed to be 
ideal up to saturation flux density of 1.6T. Regardless of 
efflorescence flux in the air gap. The force exerted on the 
moving parts is computed as a function of the air gap y. Next, 
effect of the dynamic equations is also checked. 
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Introduction 

Two key factors in selecting an actuation mechanism 
are reliability and cost. In terms of reliability, 
eelectromagnetic actuators can be controlled with low- 
cost electronics. In terms of reliability, relay design 
requires the creation of executable models for both 
features, the elastic spring and the tension generated 
by electromagnets during the performance measure 
the relative energies associated with both, 

displacement coil and magnetic attraction for each 
component. Than another, provides a range of design. 
Solving nonlinear ordinary differential equations with 
two couples Accurately simulates the dynamic 
response characteristics of the instrument. Here, less is 
studied to the elastic spring. And more, the 
electromagnetic equations, force calculation, and then 
examines the dynamic equations. 

The studies done on this field including: (1) Inductance 
of a plunger-type magnet: In this paper, a method for 


the evaluation of this function, using Mathcad, is 
proposed and two verifications of the formula are 
performed: using a finite element method and 
experimental determinations of the total inductance for 
large range of the gap variation. The experimental 
values are compared with the sum of internal 
inductance, calculated with formula and external 
inductance calculated with finite element method 
[Cividjian et al. (1998)]. (2) In second paper: Equations 
derived from the continuum design sensitivity analysis 
(CDSA), in conjunction with the material derivatives 
for a continuous medium and using the energy-based 
approach, have been successfully applied to the 
calculation of both total force and force distributions. 
The resultant expressions are similar to the Maxwell 
Stress Tensor, Magnetic Charge method, and Virtual 
Work Method but have several advantages over the 
traditional approaches. Numerical implementation of 
the scheme leads to efficient calculations and improved 
accuracy [Kim et al. (2007)]. (3) The third paper: It 
present a novel moving-iron, high-pressure, high- 
speed electromagnetic actuator that utilizes permanent 
-magnet (PM) shielding for increasing air-gap flux. We 
analyze its static and dynamic characteristics by using 
the finite-element method, taking into account the 
nonlinear characteristics and the eddy-current loss of 
the magnetic material. The experimental and simulated 
results agree well and show that the actuator has a 
displacement of 0.6 mm, closing time of 2.24 ms, and 
opening time of 7.78 ms without latching force [Man et 
al. (2010)]. Therefore, in this paper, tried that. 
Precautions be followed. 

Calculations of the Electromagnetic Force 

The mechanical attractive force generated by passing 
current through a solenoid can be determined using 
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Ampere and Faraday's laws of magnetism. Ampere's 
law stipulates that the integration of magnetic field 
intensity, H, around a closed contour of length 1 is 
equal to the net current crossing the surface of the 
closed contour: 

j)H.dl=Ni (1) 

where i is the applied current and N is the number of 
times the current encircles the contour. The integral of 
the field intensity in the contour is equal to the total 
magnetic flux, cp , divided by the cross sectional area. A, 


and the magnetic permeability, A , of the material: 


H = 


/iA 


= Ni 

* /jA 


( 2 ) 


(3) 


where the flux can be separated from the integral by 
defining a term for the magnetic resistance called the 
reluctance, R. Ampere's law can then be rewritten in 
terms of the magnetic flux and resistance through a 
material as: 

Rep = Ni (4) 


The reluctance of air gap: 

R = Xgap ~ X (5) 

/u A 

The reluctance of the air gap with permeability, u 0, is 
defined by the displacement of the spring, x, from its 
initial gap distance, x g a P . Ampere's law was developed 
for time invariant fields. Thus when considering only 
Ampere's law, x=0 and the reluctance remains constant. 


TABLE I PHYSICAL CONSTANTS FOR ELECTROMAGNET DESIGN 


Physical constants 

value 

Applied Current (A) 

3 

Number of Turns (N) 

too 

Permeability of Free Space ( Ao ) 

4ttx1 O^N / A 1 

Permeability of Nickel ( A ) 

800 x ju fj 

Area cross plunger(m2) 

3.8*10 -> 

Initial air gap(mm) 

0.5 


The flux linkage is defined as the number of turns in 
the coil multiplied by the magnetic flux passing 
through the cross-sectional area of the coil. The voltage 
may then be derived as: 

y _ d ( Ncp ) _ d(Ncp) di + d(Nep) dx 
dt di dt dx dt 

Equation (6) can be rewritten using the inductance, L, 
of the circuit times the current, and the speedance, Ks, 
of the circuit times the velocity of the plunger: 

V=L(x,i).—+K(x,i).— (7) 

y ' dt s dt 

The speedance is the coefficient of motional EMF 


across the coil and can be derived as: 


K,(x,i) = N^L = N l- 
dx dx 


Ni_ 


( 8 ) 


Assuming no electrical power loss in the coil, the 
change in work done as the total flux changes is: 

W = J i.d(Ntp) = j Rcp.dcp = ^iNtp 


(9) 


F = 


( 10 ) 


The force acting on the second magnetic element is: 
dW _ Ni dtp 
dx 2 dx 

Substituting for the reluctance in the solution of the 
magnetic flux, the force provided by the electromagnet 
on the magnetic plunger can be written as: 

(M) 2 1 


Fmag 2 u 0 A 


x g °e~ x 

Fo A 


( 11 ) 


The Exact Dimensions of the Model 

Dimensions of the simulation model in the software 
Maxwell Presented as follows: 



TABLE II DIMENSIONS OF THE SIMULATED MODEL 


Var. 

mm 

Var. 

mm 

Var. 

Mm 

a 

40 

u 

2.5 

R 

10 

c 

9.5 

V 

2.5 

a' 

40 

b 

10.5 

Lt 

9 

hb 

39 

di 

0.5 

3 

0.5 

hi 

10 

D2 

0.5 

s. 

0.5 

h2 

10 

ch 

0.5 

h.4 

60.5 

t 

4.8 

h 

60 

Li 

10 



L 

60 

L,2 

10 




Analysis and Simulation 

In this section plunger simulated by two state. 
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Sinusoidal Voltage 


For example, we have applied sinusoidal voltage 
V(t) = lOCh/2 cos(120;rf) . displacement to 1mm we 
calculate amount of force applied to the piston. 
Current is: 


i{t) = 


VO) 


100V2 cosl20 nt 
jr 2 +(Lcof 


T n 2 u,,am 100 2 x4^x10 7 x380x10 6 x120^ 

Leo = z = r 

2x 1x10 3 

.. . 100^2 cosl20;rt IO0V2 

i(t) = , = cosl20;rt 

V4 2 +(1.8) 2 4.39 

1 400V2 cosl20;zt., 

/ = -( )"x... 

2 4.39 

100 2 x4;rxl0~ 7 x3.8xl0^ f 1 'l 

2 t 10' 6 J 


/ = -1238 -89cos 2 120^ = -1238-89 


l + cos240/rt 


=619 -4451V 


In the figure 2, Displacement one mm moving parts is 
shown. 



FIG. 2 DISPLACEMENT ONE MM MOVING PARTS 



FIG. 3 ELECTROMAGNETIC FORCE, FOR A SINUSOIDAL 
VOLTAGE V(t) = 1 0()V2 cosO 2()7it) 


In the figure 3, the peak power is 532.8N with the 
calculated value 619.44N is about 14% difference. This 
large variation due to the reluctance of the fixed and 
removable nozzle gap is negligible. So it is better to be 
considered. Positive sign, because it is also the size of 
the force, the plot here. 


In the figure 4, peak current is In the shape about 
34.27A. With calculated amount 32.21A, a difference of 
about 6%. 



DC Voltage 


For a voltage Vdc=12V, and the flux density of 1.6T 
assuming that the magnetic material is ideal. To 
calculate the amount of force to Displacement 1mm, 
Current steady state is as follows: 

. V 12 ,, 

1 = — = — = 3 A 
r 4 

cp = BA = 1.6x3. 8x1 0" 4 =6.08xl(T 4 vtT> 

_ _100x3 
<P 




<P 


, N 

, = ,0439xl0 6 — » Z = = 0.228 H 

6.08x10^ R 


y = R/.t 0 A = .0439 x 10 6 x An x 10“ 7 x 3.8 x 10“ 4 




= 2.1x10 ‘ =0.021 mm 

100 2 x4;rxl0- 7 x3.8xl0- 


10- 


= -21.488 N 


This means that that for every y <.021 mm, the force 
exerted on the piston due to the saturation of magnetic 
material, is fixed. The gap is greater, R is greater 
reluctance, L is less inductance, B is less; then magnetic 
material is not saturated. And vice versa. 


Figure 5 compares the calculated value of the force 
difference is the difference between taking the spring 
constant K occur in software. 


We see in Figure 6, the peak current is 3A. And is 
accordance with calculated value. 
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Winding 


Maxwell2D Design 1 



FIG. 6 THE RED COLOR REPRESENTS THE CURRENT 


Dynamic Equations 


Using equations previously described. Coupled set of 
first order nonlinear differential equations can be 
developed to make a completely model tool. 


dx 

dt 


m.— = F g —K.x—b(x,i).v+K s (x,i).i 


( 12 ) 


L{x,i).— = V - K (x,i).v - R..i 
dt 

Which, b is the damping coefficient of viscosity, and K 
is mechanical spring constant. 


, 3 .v'.A 2 

b = r 

2 .n( x gap -x) 

Which, v' — 1.583xl0 -5 kg / ijn.s') 
b = 2.872N.S / m . 


( 13 ) 

A = 3.8xl(T 4 «r then 


Inserting values with K s - .2 , K = 5335 , L = .228 , V a = 

1 2vDC , in = .009/ty , F g = 0.09/V and writing equations 

in the form of state equations and solve it in software 
obtains the required waveforms. 

i, = x 2 


x, =— (0.09-5335.X. -2.87.x 2 + 0.2.x,) ^ 
~ .009 

x, = — - — (12 -0.2.x 2 -4.x,) 

0.228 - 


These three equations are solved in MATLAB using the 
solver ODE45. 

In the figure 7, after connecting the power piston may 
fluctuate. And in the end proved to be the .1 mm. But 
the position is not zero. 

In the figure 8, the piston is fixed and have no speed. 

In the figure 9, Current, according to calculations at 
steady state, is fixed to the 3A. 

But, with a alternating source Position changed 
oscillatory. Due to the structure Plunger, here is 
unacceptable. But it can be use in a tool as a Moser 
machine that is as small fluctuations, but the two sides 
are allowed to fluctuations reciprocating this AC 
voltage is applied. Note to Figures 8 to 10. Of course, in 
theory, be made condition on the position that never 
becomes negative but this condition also affects the 
other two state variables. 

In the figure 11, at the moment, connect the source; in 
piston speed, several peak occurs. This is not too far 
from expectation. 

The coefficients of the equations were estimated. To 
achieve more accurate response parameter should be 
calculated as the exact Ks. Formalin calculate 
parameter such as b, also are approximate. Other state 
variables like the electromagnetic force, energy field or 
flux linkage can be added into the state equation. But 
we could directly use Rung-kuta format to solve the 
equations. Again, this method is based on ODE45. 



FIG. 7 PISTON POSITION CHANGES 



FIG. 8 PISTON SPEED CHANGES 



FIG. 9 CURRENT IN THE COIL 
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FIG. 10 PISTON POSITION CHANGES 




FIG. 12 CURRENT IN THE COIL 


Conclusions 

Seen, the calculated values and simulation values were 
somewhat close together. If, in the calculation, 
regardless of the approximation, the accuracy is higher. 
But the possibility of lack of convergence, also not 
unexpected the correct answer in Maxwell Software. In 
the analysis of the dynamical equations we assumed 
constant the L and R. Ideally, of course, will mean that 
B=1.6T. But in the General state changes of R and L, 
also changes the equation. The current form will have 
minor changes in the transient state. Dependence of 


dynamic and the electromagnetic equations causes to 

solve the range of problem and design. 
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